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Abstract. As starburst galaxies show a star forma- 
tion rate up to several hundred times larger than the 
one in a typical galaxy, the expected supernova rate 
is higher than average. This in turn implies a high 
rate of long gamma ray bursts (GRBs), which are 
extreme supernova events. We present a catalog of 
127 local starburst galaxies with redshifts of z < 0.03. 
Using this catalog we investigate the possibiUty of 
detecting neutrinos from Gamma Ray Bursts from 
nearby starburst galaxies. We show that the rate of 
long GRBs is correlated to the supernova rate which 
in turn is correlated to the far infrared output. For 
the entire catalog, 0.03 GRB per year are expected to 
occur. The true number can even be higher since only 
the brightest sources were included in the catalog. 
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I. Introduction 

In our work fTl we present a systematic investigation 
of radiative emission from starburst galaxies. We ex- 
plain there the correlation between radio emission and 
emission in the far infrared and high energy photons 
and neutrinos. In this contribution we will focus on the 
possibility to detect secondaries from gamma ray bursts 
in starburst galaxies. In section the sample of 127 
starburst galaxies used will be presented. Section |in] 
explains the connection between starburst galaxies and 
gamma ray bursts while section |lV]presents the neutrino 
spectrum for GRBs. Expected event rates in the IceCube 
detector for GRBs from nearby starburst galaxies are 
shown in sub section IIV-AI 

II. A LOCAL STARBURST SAMPLE 

First a sample of 127 nearby starburst galaxies is 
seleected. For the data of the individual sources see 
H] Appendix A. This is a local sample containing only 
sources with z < 0.03. The sources were selected from a 
larger catalog of starbursts which contained 309 sources. 
To ensure completeness of the sample, cuts in the radio 
flux as well as in the FIR flux were applied. To ensure 
that the selected sources were indeed starbursts, only 
sources with a high ratio of FIR flux and radio flux, 
<S'60ai/'5'i.4GHz > 30 were selected. Sqo^ denotes the 
FIR flux density at 60 fj,m wavelength measured by the 
IRAS satellite, S'i.4ghz denotes the radio flux density 
at 1.4 GHz. This criterium removes possible contami- 
nation of the sample by Seyfert galaxies. The ratio is 
~ 10 for Seyfert galaxies and ^ 300 for starbursts. 



Further, luminosity cuts were applied. It was required 
that S'i.4GHz > 20mJ}Q and Sqo^ > 4Jy. 

Figure [Ta] and figure [lb] show the luminosity of 
the starburst galaxies versus their luminosity distance 
at 60 /im resp. at 1.4 GHz. The crosses represent all 
309 starbursts that were selected in the beginning, the 
triangles show the remaining 127 sources after the cuts 
Seofj > 4Jy and S'i.4ghz > 20mJy and z < 0.03. 
The dotted lines represent the sensitivity for 4Jy and 
20 mJy, respectively with those cuts applied, a complete, 
local sample in both wavelengths, FIR and radio is 
obtained. Since all sources in the sample have a redshift 
of z < 0.03, many of them are located in the super 
galactic plane. Therefore, the number of sources with 
a flux density larger than S, N{> S) is expected to 
follow a powerlaw behavior of S^^ — 5^^^. For a 
flat cylinder a pure S^^ behavior is expected while 
a spherical distribution results in an S^^ '^ behavior 
Figure |2] shows the logarithmic number of sources above 
an FIR flux density Sqq^. The data was fitted with the 
following function: 



iV(> S)^No-{S + Soy^ 



(1) 



Here, A'o, 5*0 and /3 are the fit parameters. Using an error 
of Vn, the parameters are determined to 

No = 3155 ± 1297.9 
So = (10.56 ± 3.78) Jy 
P = 1.2±0.2. 

The behavior iV(> S) ^ S*^^-^^*^'^ matches the expec- 
tation that the function should lie between 5^^ ° and 

The next section will explain the connection between 
starburst galaxies and gamma ray bursts which will result 
in a flux estimation for high energy neutrinos from GRBs 
inside starburst galaxies. 

III. Gamma Ray Bursts and starbursts 

Starburst galaxies show an enhanced rate of supernova 
explosions due to their large star formation rate. Thus 
an increasing rate of long Gamma Ray Bursts (GRBs) 
directly linked to SN-Ic events [2] is expected. If long 
GRBs are the dominant sources of UHECRs, the contri- 
bution from nearby objects should follow the distribution 
of starburst galaxies. In the following calculations, it is 
assumed that every SN-Ic explosion is accompanied by 
a particle jet along the former star's rotation axis, i.e. by 
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(a) Luminosity at 60 /im versu.s tlie luminosity distance 
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(b) Luminosity at 1.4 GHz versus the luminosity distance 



Fig. 1. Luminosity-distance diagrams. Crosses denote all 309 pre-selected starburst galaxies, triangles show the remaining after the cuts 
SeOfi, > 4Jy and S'i.4GHz > 20mjy and z < 0.03. The dashed lines show the sensitivity for Seo^ = 4Jy and 51.4 0112 = 20mjy 
respectively. 
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Fig. 2. log{N) — log{S) representation of the catalog. An 5 ^'^-fit 
matches the data nicely, with a turnover at 5o = 10.56 Jy. 



a GRB. The opening angle of the GRB jet, 9 determines, 
how many SN-Ic can be observed as GRBs, see e.g. [3!, 



11. 



"GRB = e • «SN-Ic 



(2) 



Here, noRB is the GRB rate in a galaxy and e = 
(1 — cos(6')) is the fraction of SN-Ic which can be seen 
as GRB. The opening angle is expected to be less than 
^ 10° for the prompt emission. Afterglow emissions and 
precursors can have a larger opening angles |5|. Putting 
the focus on prompt emission, an optimistic opening 
angle of 10° is used, yielding 



e = 0.015. 



(3) 



Further, observational data show that core collapse su- 
pernovae of type Ic contribute with 11% to the total 
supernovarate in starbursts |l6|. Using equation |2] the 
GRB rate in a starburst galaxy is directly correlated to 
the supernovarate fisN, 



«GRB = e ■ C ■ "SN 



(4) 



with ^ ~ 0.11 as the fraction of heavy SN among all SN. 
The supernova rate is correlated with the FIR luminosity 



of the galaxy ||7l, 

nsN = (2.4 ±0.1) • 10 
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(5) 



The FIR luminosity is expressed in terms of the solar 
luminosity Lq = 3.839 • lO'^'^ergs and is given in the 
range of 60 fim and 100 /im by 181 



(6) 



Here, di is the luminosity distance of the individual 
source and 

Ffir = 1.26- 10-14. (2.58 . S^o/^ + ^loop) Wm^ (7) 

is the FIR flux density as defined in |9|. Sqq^ and 5ioo/j 
are the measured flux densities at 60 /im and 100 /im, 
both measured in Jy. Using equation |5] to determine the 
supernova rate, equation |4] yields a rate of 



"-GRB = 3.8. 10 
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(8) 

per starburst. For a 1 km"^ neutrino detector with a life- 
time of 10 years (like IceCube or KM3NeT) luminosities 
of around 3 . 10^'^ . Lq ^ 10"*^ ergs are required for a 
single event within these 10 years. None of the sources in 
the catalog provides such a high luminosity. However, if 
a larger number of starburst is considered for an analysis, 
the total luminosity increases and so does the probability 
of observing a GRB. Figure |3] shows the total GRB rate 
for a number of A^starbursts galaxies. 
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(9) 

In the figure the GRB rates achieved in the single star- 
bursts are summed up, starting with the most luminous 
source, adding sources in descending luminosity order 
The points show the GRB rate summing up over all 
starbursts in the sample, starting with the strongest one 
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discussed in detail in |[T]. Their numerical values were 
determined to 



nr 



Fig. 3. Number of GRBs per year in the starburst catalog, including 
A'atarbursts sources. Starting with the strongest one. The total GRB 
rate in the sample, including all 127 sources is 0.03 yr"\ this means 
that a GRB could be observed every 30 years on average. The 
total GRB rate just in the northern hemisphere is 0.02 yr~^ or an 
occuiTence every 50 years, the data shown in squares. These source 
lie in IceCube's FoV. 



and adding the next strongest sources subsequently. On 
total, 0.03 GRBs per year are expected to be observable 
in the sample. The squares display the total GRB rate, 
summing up sources in the northern hemisphere, which 
corresponds to IceCube's Field of View (FoV). Here, 
0.02 GRBs per year are expected. This number can be 
enhanced significantly when taking weaker sources into 
account which were not included in the sample in order 
to ensure completeness, see section HIl 

IV. Enhanced neutrino flux from GRBs in 

STARBURSTS 

Due to the high atmospheric background seen by high- 
energy neutrino telescopes, the detection of a diffuse 
neutrino signal from GRBs in nearby starbursts will not 
be possible. However, with a timing analysis one might 
be able to identify Gamma Ray Bursts in neutrinos. In 
such an analysis, the location of a nearby starburst can 
be chosen as a potential neutrino hot-spot. By selecting 
a time window of the typical duration of a long GRB 
(~ 100 s) the atmospheric background can then be 
reduced to close to zero. In this context the general 
neutrino intensity and in particular the possibility of 
neutrino detection with IceCube are discussed. 
For the first time the neutrino energy spectrum during the 
prompt photon emission phase in a GRB was determined 
by Waxman&Bahcall ifTOl . ifTTI and can be expressed as 




< E < 



(10) 

The spectrum includes two spectral indices, and /3y, 
two break energies, e\, and ef and a normalization factor 
A,^. For the GRBs in starbursts, these parameters were 



A. 
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The normalization constant is calculated for each 
individual source. It depends on d^^ as well as on 
the fraction of energy transferred into electrons and the 
fraction of energy transferred into charged pions. In ad- 
dition, the normalization of the neutrino spectrum scales 
with the luminosity of the burst. This released energy 
varies from burst to burst. In addition to this burst-to- 
burst fluctuation, regular GRBs are distinguished from 
low-luminosity bursts. Regular, long bursts emit a total 
isotropic energy of 10''^ erg for a duration (fgo) of the 
burst of w 10 s. Low-luminosity bursts last longer and 
and have a lower luminosity. Although only few low- 
luminosity bursts are observed yet, they are expected 
to be much more frequent than regular GRBs. For this 
class, we expect an energy release of ~ 10^*^ erg within 
around 1000 s. The closest burst observed so far was 
GRB980425, which was found to be associated with the 
supernova SN1998bw [fT2l . The host galaxy lies at a 
redshift of only z = 0.0085. This burst shows a total 
energy release of ^ 10^^ erg, which is an extremely 
low-luminosity burst. As the luminosity distribution is 
not well-known at this point, due to low statistics, we 
use a fixed value of 10^^ erg. An actual burst can be 
about one order of magnitude more or less luminous. 
Now, to estimate the neutrino flux from a standard GRB 
for a single starburst in the sample, dependent on the 
distance of the starburst, the normalization is calculated. 
The other parameters are kept constant and hence the 
results can only serve as a rough estimate. Both, the 
break energies as well as the spectral indices vary for 
each individual burst as described in 111 31 . 

A. Expected event rates in IceCube 

As shown in figure [3] 0.02 GRBs per year are ex- 
pected to occur in the 96 starbursts in the sample in 
the northern hemisphere, IceCube's FoV. However, this 
rate can be enhanced if all sources in the northern hemi- 
sphere would be taken into account. For completeness 
reasons only the brightest ones were considered here. 
This enhances the possibility to detect a GRB from a 
starburst in the super galactic plane within the lifetime 
(10 years) of IceCube. The prospects for KMBNeT are 
slightly worse, since only 0.01 GRB per year is expected 
in the southern hemisphere, where KM3NeT's FoV will 
be focused. 

The number of events per GRB expected in IceCube 
can be calculated by folding IceCube's effective area 
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Fig. 4. Histogram of tlie number of events in IceCube from the 
96 starburst galaxies in the northern hemisphere. Depending on the 
distance of the starburst, a burst would result in between 0.1 and 
several 100 events in IceCube in a small time window of 10 — 100 
seconds. A regular burst at ^ ~ 1 — 2 gives only 0.027 events as 
indicated in the figure. The main reason for the increased signal is 
that the bursts would come from starbursts closer than z = 0.03. 



Aes El with the GRB spectrum dNJdE^ 



(11) 



Here, the weak dependence of the A^ff from the dec- 
lination of the burst is neglected. For the threshold 
energy i5th — 100 GeV is used. This is the general 
detection threshold of IceCube |15|. Since events can 
be selected in a small time window with the typical 
duration of a long GRB, 10 — 100 s, the atmospheric 
background can be reduced close to zero. Figure |4] shows 
the histogram of the numbers of events expected in 
IceCube from an average GRB with an isotropic energy 
of E}^° = 10^^ erg from the 96 starbursts in the sample. 
These numbers range from 0.1 to about 300 events per 
burst, depending on the distance of its host galaxy. These 
numbers lie between 1 and 5 orders of magnitude above 
the numbers of events for GRBs typically observed by 
satellite experiments like Swift, BATSE and Fermi. If 
such GRB occurs in our sample of starburst galaxies it 
is a unique opportunity to study the hadronic component 
of GRBs. 

V. Conclusions and Outlook 

A sample of 127 nearby starburst galaxies has been 
presented, where we selected sources with z < 0.03 
and radio (FIR) luminosities above S'eo^ > 4Jy and 
5*1.4 GHz > 20mJy. 



The expected event rates for high energy neutrinos in 
the IceCube detector have been calculated to reach up to 
more than 100 events per year. If a GRB occurs in one 
of the 127 starbursts, for detectors with a wide field of 
view the detection probability is best. IceCube should 
be able to see the strongest events 100 neutrinos) 
immediatly, while the weaker events {^-^ 5—10 neutrinos) 
can be found in a specific analysis for the selected 
starbursts. Since most of the sources are located in the 
northern hemisphere, detectors like IceCube, HAWC and 
Auger north are optimal for such a study. A detection 
of neutrinos from a GRB would enable detailed studies 
of hadronic emission processes of GRBs. 
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